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1. INTRODUCTION 


In various branches of applied mathematics one is often confronted 
with problems which involve a certain function of the time, called 
the resultant W (t) in what follows, that can be looked upon as built up 
through a cumulative process started and continually kept up by 
another function of time, the originator n (t), the effect of which on 
W (t) is brought about through the intervention of a certain function 
of elapsed time, the distributor w (€). 

The following are some instances drawn from the field of population 
statistics: The total existing population, or certain subgroups of the 
population, such as the number of men, or women, or men in certain 
age groups, at a certain moment of time, or the annual number of 
deaths or of marriages, are all resultants, i.e., functions of time which 
are originated through another function of time, namely the number of 
births in previous years; in these examples, the distributors are certain 
functions derived from the age distributions at death, at marriage, etc. 
Again, the total reserve at a given moment for a particular kind of 
insurance is a resultant which is originated through the “production,” 
i.e., the new insurance written each year; here the distributor is the 


1 Received June 8, 1928. The present paper is a revised and somewhat generalized 
version of the author’s article Ueber periodische und angendherte Beharrungszustdnde, 
Skandinavisk Aktuarie Tidskrift. 1926. 
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function of elapsed time that represents the total reserve resulting from 
a unit of production of one year for the kind of insurance considered.? 

In all the above cases the effect of the distributor is limited by the 
maximum length of human life. An analogous limitation is present in 
a great variety of cases occurring in other fields. 

The most important case is perhaps that in which the cumulative 
process is such that the originator-element n (r)dr for the moment of 
time t = 7 is the cause of a contribution of amount w ()n (r)dr to the 
resultant at the subsequent moment of time t = 7 + & We then 
have, since w (t) = w(t — 7), 


t 
W (é) = i) w (t — r) n(r) dr SRR: PN a Ree acrn (1.1) 
As, however, the originating process will in general start at some 


definite time which we may select for t = 0, we may assume n (r) = 0 
for r < 0; thus we obtain 


t 
W (2d) = ('wd-Dn@as iis anlar ak ctochins anal’ (1.2) 
If the originator n (t) is a constant n, we have 
t 
Wi) =n wll — Adres. (1.3) 
which, upon substituting  — +r = §£,dr = —dé (since ¢ is a definite 
fixed time), may also be written 
t 
W (2) =n{w@de Veal ee ou, 2 Sk (1.4) 
If the effect of the distributor is limited, i.e., w (¢) = O for é > T, 
where 7’ is some positive constant, we get for values of t > 7’ 
T 
W (i) = nf w (g) d& = const.............. (1.5) 


This last equation expresses a well known equilibrium theorem. 


* The advantages of considering actuarial reserves from this point of view has been 
pointed out by the writer in an article in Skandinavisk Aktuarie Tidskrift, 1921, p. 210. 
See also the work of Norwegian actuaries in Nyt beregningsgrundlag for livsforsikring, 
Oslo, 1922, p. 27, and Fr. Lange-Nielsen, A proposal for a new basis . . . . , Skand. 
Akt. T., 1922, p. 235. 























In practical cases where the originator and the distributor are 
known with sufficient accuracy, these equations give a means of 
forecasting the magnitude of the resultant. In particular, if n = 
const., and the distributor is limited, equation (1.5) gives a means of 
predicting the time from which there will be equilibrium in the sense 
of W = const. It is obvious that it would be of great practical impor- 
tance to obtain criteria for cases in which equilibrium will be approxi- 
mately reached when the originator n (¢) is not constant; the attempt to 
do this is the chief aim of the present paper. 


2. THe ORIGINATOR AND THE DISTRIBUTOR AS SOLUTIONS OF 
INTEGRAL EQUATIONS 


In the following we shall write equation (1.2) in the form 


W (i) = [went —tdt = [we — 1) n(r) dr...(2.1) 


A generalization of this equation would be 


W (i) = } z wi; (€) ni (¢ — &) dé.... 2.0... (2.2) 


where w; and n; are two series of distributors and originators; if either 
the w, or the n; are identically proportional, this case reduces to (2.1). 
In the present paper, attention will be confined to the case (2.1). 
Equations (2.1) are integral equations, in w () and in n (r) respec- 
tively, of the so-called generalized Abel type; to obtain practical solu- 
tions of them, it is desirable (when the nucleus remains finite) to 
transform them by differentiation with respect to t, whence we obtain*® 


[wo@nt-ode 








W' (t) _ 

(0) = w (t) + r= 0) Eb ote ele ele (2.3) 
W’' (t) _ [n@wt-sas 

ster" Oi 7) aa (2.4) 


If in (2.3) the quantity n (0) vanishes, the differentiation must be 
repeated until n” (0) + 0; and similarly if w (0) vanishes in (2.4). 
Equations (2.3) and (2.4) are integral equations, in w and n, respec- 


3 E. B. Witson, Advanced Calculus, p. 283. 
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tively, of the so-called Poisson type; the existence theorems established 
by Volterra‘ give a solution of the form 


we [re-ow wae 


w(t) =) n 0) 





or, transforming the variable of integration from ¢ to r = ¢t — &, 
dr = —dé, 
t 


(he ; 
a n (0) 





in which Tf (¢ — &) is called the resolvent kernel or solving function; 
numerical methods for determining [ have been given by Whittaker,® 
Kameda,‘ and others. If the kernel n’ (¢ — £)/n(0) = K(t — £) in 
(2.3) is of the form K (x) = 2; a; exp (civ), where the a; and c; are 
constants, I will be of the same form [ (2) = 3; a; exp (7; x), where 
the a; and y; are simple functions of the a; and c;. In particular, 
if n (t) = a + 6 sin d t, K (t —£) = 2b cos d (t — &), where a, 8, 
r, 6, are constants, we have 


Tr (t) = 2 be” (cosy it — vsiny 2d), 2 > Be 
r (t) = be“ [(1 — v) e** + (1 +) 7"), x2 < Be}... (2.7) 
p= by = + Vp? — BF 


3. Some GENERAL RELATIONS BETWEEN THE FLUCTUATIONS OF THE 
ORIGINATOR AND THE VALUE OF THE RESULTANT 
AT A GIVEN INSTANT 


Let W, (¢) be the function which expresses what the resultant would 
have been at time ¢ (with the given form of w) if the originator had 
been equal to a constant m; and let us find an expression for the 
deviation U (t) = W (t) — Wo (t). 

If m = mn (0), we have from (1.4) and (2.6), remembering that 
W (0) = O and that the value of a definite integral is independent of 
the particular variable in terms of which the integrand is expressed, 


* Goursat, Cours d’ Analyse, 3 ed., T. III, chap. xxx. 

‘WuitTaKeR. Proc. Roy. Soc. Lond. 94: 367. 1918. WuxiTTakerR and Rosinson, 
The Calculus of Observations, p. 376. 

* KamepA, Tohoku Math. Journ. August, 1924. 
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Um =WW-(n@wwat 


= far Vrowe-oa 


or, interchanging the order of the integrations by Dirichlet’s Formula,’ 


uw =larfr@wG-da- 


Hence, evaluating the second integral (in which W’ is now the deriva- 
tive with respect to 7), we have 


U@=(r@we-Hde=l'ru-pw wa 
Again, from (2.1) we get directly 
Ui) = ['we@ne-Hde—m (weds 


The deviation U (t) may be computed from either of these equations 
for any value of ¢t. The last equation is also valid if n is an arbitrary 
constant. 

Generally, let W (t) and H (t) be two resultants, generated by the 
originators n (¢) and A (t) respectively, and with the same distributor 
w (€). The function A (t) will be called the trend function. The 
difference u (t) = n (t) — h (t) will represent the deviations or fluctua- 
tions of the originator n (t) measured from the trend function, and 
U (t) = W () — Ai (t) will represent the corresponding deviations of 
the resultant. 

We shall first suppose h (t) to be an arbitrary function without any 
connection with n (t). By (2.1) 


vw _ jut-du@ae 
aliens frd-dw@a: 
0 


UIE OR yoo 
liters h(t — &) w (&) dé 


[na-ow@de 


7 Goursat, Cours d’ Analyse, 4ed., T. I., p. 299. 
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Thus, by a well known property of definite integrals, 











tlu (t — €) 
2.) Slen=alh- 2° ele 
H(t) “ 
J rd-ow@ae| 
Hence 
UW) [i|ac-owe@lae 





<M 
H (t) Lf 
[Jee - owe ae| 


where M is the maximum of | u (r) /h (r) | in the interval 0 < r < 
If h (t—£) w (€) has the same sign throughout the interval 0 < & < 
which in practice it usually will have,* we have simply 


U (t) 
Hi) <M 





Therefore we have ‘ 

TueoreM I. Whatever the trend function h (7), if the product 
h (r) w (t—7) does not change sign in the interval 0 < - < ¢, the 
relative deviation of the resultant W (¢) from the value which this func- 
tion would have assumed if the originator had been the trend function 
h (r) has at any given point ¢ an absolute value not greater than the 
maximum absolute value of the relative deviation of the originator 
from the trend function in the interval 0 < 7 < t. 

In other words, since this theorem holds for any ¢, a damping of the 
relative fluctuations of W (t) as compared with the relative fluctua- 
tions of n (¢) must in all such cases take place. Upon the intensity of 
this damping will depend the possibility in practice of computing 
W (t) as if n (t) were a suitable trend function h (t) of a less compli- 
cated form than n (¢) but still representing the main features of n (t)— 
e.g., as if n (t) were a constant in cases where the underlying trend of 
n (t) is recognized to have the form of a horizontal line. 

The damping effect will be more fully examined below, through the 
introduction of the standard deviation as a measure of the fluctuations 


of W (¢). 


8 In most cases, both h (t) and u (¢) will be positive. 
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If the originator n (¢) has a marked variation confined to a short 
interval in the neighbérhood of time ¢, and if the distributor w (z) is 
positive and has a pronounced maximum at x = a, then there will 
obviously be a more or less pronounced fluctuation in the resultant 
in the neighborhood of time ¢ + a. We may obtain some information 
about this effect by considering the originator to be increased over a 
small finite interval of length p = 29, say: Let m; (x) be a function 
which is positive in the interval 0 < x < p, and equal to zero outside 
this interval; nm, (x) might be called the “extra originator.’”’” We 
wish to know the effect on the resultant W when, starting from any 
time ¢, the extra originator is added to or subtracted from n. 

The resultant which would be produced at any instant of time 
t+y >t+pbymalone, viz., 


Wi (y) = fw fal Om OER go. conan (3.1) 


might be called the ‘extra resultant.’’ The problem is to determine 
the maximum of this function. If there is a maximum at y = 7 + q, 
7 must satisfy 


+@ 
fw @-dm@tode=0 


if w (x) is symmetric about its maximum x = a, and n; (x) symmetric 
about the midpoint x = gq, this equation will be satisfied for » = a. 
Hence, the maximum effect on W of a temporary, symmetric fluctua- 
tion in the originator n is produced after a lapse of time which, counted 
from the midpoint of the fluctuation in n, is equal to a, where a is the 
argument for which the symmetric distributor w (x) has its maximum. 

The conditions under which this solution holds are often realized in 
practise, because several distributors encountered in statistical and 
actuarial work, e.g., are approximately of the symmetric type. A 
striking instance is the reserve 


( ) L, +t 
w=c\l— bas Sakis ] 
a, z 


of a sufficiently great number of whole life insurances entered at age zx. 





4. Tae Dampine CogEFFICIENT 


A formula will now be established that expresses the way in which 
the general character of the fluctuations of W (t)inz —q@ <t<z+4¢q 
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depends on the fluctuations of n (¢) in the time previous to z + q. 
The distributor will be assumed limited in its effect, i.e., w (x) = 0 
forz > T; hence fort > T 


T T 
W (it) -{ winit—ode= fi w (t — ¢) n (&) dé.. (4.1) 


Assume also that z > 7 + q; T is the instant of time at which equi- 
librium would have been reached, i.e., the time from which W would 
have been constant, if the originator n had been independent of time, 
and it will be called the equilibrium point regardless of whether or not 
equilibrium actually is reached there. The equilibrium point evi- 
dently is independent of the constant value assumed for n, and depends 
only on the distributor w (x). 
The function 


w (x) 


w (x) = aR 
J, w@ae 


will be called the normalized distributor, or the distribution function. 
We further introduce the means 


1 (zt+¢ 1 (zt¢ 
h (2) -7 J n(r)dr; H(z) = a ae W (r) dr 





in which p = 2q; if z is varying, h and H can be regarded as the moving 
averages of nand W.* The deviations of n and W from their means 
will be denoted by 


n (t) — h (2) 
W (t) — H (2) 


u (t, z) 


U (t, 2) 


It should be noted that the position of the interval which determines 
the means from which u and U are measured is fixed independently 
through z, and does not change with ?. 


* In my original analysis of the character of the oscillations of W (¢) over a given 
interval, I considered only the case of a periodic originator, which evidently amounted 
to considering the trend function A to be a constant. The more general definition of 
hand H, above, as moving averages has been suggested to be by my friend and colleague 
Dr. Frisch. 
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We evidently have for any z 


z+q z+q 
i) w(r,2)dr = f U (r,z)dr =0 


—— 


The means and the deviations satisfy the same equations as do the 
func:ions themselves; in fact, we have by (4.1) 


T 
H(z) -{ w (t) h (2 — ) dé............. (4.2) 
and 
T 
U (t, 2) -{ w (€) u(t — &, 2 —&) dé.......... (4.8) 
fort > Tandz>T7 +4. 


The standard deviations s (z) and S (z) of n (¢) and W (t) over the 
intervalz —q <t < 2 + qare defined by 


1 z+¢@ 
#() => f u? (r, 2) dr 
Ddz-a 


1 (zt+¢ 
S? (2) = i) it a 
° P es 


Finally, we introduce the coefficient of correlation C (25, z) between 
the 25-spaced ordinates n (¢—6) and n (¢ + 64) of the originator, the 
correlation being taken over the interval z — q <t <z2+4@: 


1 z+@ 
~{ u(r +6,24+ 5) u(r —6,2—58) dr 
C(%:2) - 





s (z + 4) 8 (z — 8) ‘Aaa 
Since C is a correlation coefficient, we always have —1 < C < 1. 
Further C (0, z) = 1 for any z, which simply means that the correlation 
is perfect between infinitesimally spaced ordihates. C is symmetric 
in 6. 

The relative standard deviation r (z) = s (z)/h (z) measures the 
closeness with which the originator n (¢) fluctuates about its mean in 
the interval z + q; and the ratio R (z) = S (z) / H (z) has the same 
significance for the resultant W (t). The ratio R (z) may therefore 
be taken as a measure of the degree to which equilibrium is approxi- 
mately realized over the interval z + q. 
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Making use of (4.3), we have for R (z), after some reductions,’ the 
expression 


G*(z) 


T (-T 
Be) = Es I. I w@wladale—ale—n)O(e—2, eae dy. .(4.5) 





where 





T T 
fw@ae(ae-oae 
G(z) =< F ee (4.6) 
Tl w@he-oae 


17 
k@ =F h(z-—édé 


We can further reduce (4.5) by means of the following formula 
due to Frisch:" Let f(a,..,2,) or, for brevity, Jo be a function of 


the n variables z,,..,2%», and let j denote {.. f dz,....d%s. 


Further let m; = ie f. be the mean of f, and let 


10 The reductions consist in writing 


T T 
rina= fal d nw (E) w (n) u (r — &, 2 — &) u(r — 9, 2 — 9) 
0 0 


then transforming the resulting triple integral 


z+¢ T T 
j arf ae Poe Sam 
z-@ 0 0 

T T z+¢ 
i) arf anf Att ee 
0 0 z—-@ 


té+y 
2 


into 


and finally putting 


r=r+ 


1 RaGNAR Friscu. On approximation to a certain type of integrals. Skandinavisk 
Aktuarie Tidskrift. 1928. The formula there given is a generalization of a formula 
which Frisch had previously given; for other generalizations of Frisch’s formula see 
Steffensen. On the sum or integral of the product of two functions. ibid., p. 44. 1927. 
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of = + Vara G-™, 


of i 
pf = — 
my 



















(m; assumed + 0) be the standard deviation and the relative standard 
deviation, respectively, of f. Finally, let 


@ — mz) (Gz — M4) 





Tt? = 
Of Og (b — a) 


be the coefficient of correlation between the two functions f, and gz. 
Then we have 


(6 — a) [s 9: = fs ‘ fo. (1 + pr po Tw) 


the means m; and m, being assumed + 0. Substituting for p;, p, and 
Tze, we see that this is an identity. 
Applying this to the two functions of two variables 


f Gn) = 8 (2 — & 8 (2 — 9) 


w (@w(n) C(¢ — 9,2 - £2") 


ll 





g (é, n) 


in (4.5), and also to the two functions of one variable w () and h (z—£) 
in (4.6), we have 


7 
f s(z—£é)dé 
R (z) = D (2) “7 a athe 
f h(z—)dé 





ye< és vate 


where 


L terete (™(? p+ 
De (2) = Gee (Ooo C(e- 12-25") dean...48) 








The ratio of the integrals in (4.7) may be regarded as an average 
relative standard deviation of the originator for the interval previous 
toz +g. Wemay therefore adopt the factor D (z) in (4.7) and (4.8) 
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as the coefficient of damping, which expresses the ratio between the 
relative standard deviation of the resultant over the interval z + q 
and the average relative standard deviation of the originator over the 
interval previous to z + gq; the magnitude of D (z) measures the 
intensity of the damping. If D (z) is a small fraction, a heavy damp- 
ing has taken place. It is true that a small value of D (z) does not 
exclude the possibility that there may have taken place, within the 
interval z + q, isolated deviations of W of comparatively great magni- 
tude, but such fluctuations must necessarily have been confined to 
very short intervals of time; in practise, such extreme cases will be of 
little importance. 

The expression (4.8) enables us to recognize certain cases in which 
D (z) will be small; we shall consider, in particular, the case of a 
periodic originator. 


5. Tue Case or A Periopic ORIGINATOR 


If the moving average h is constant, i.e., if the originator is periodic 
with period p, then A, s, H, and S (for z > 7’ + q) will be constants, 
hence oy = o, = 0. Furthermore, u (¢, z) and therefore C (4, z) will be 
independent of z. The ratios r and R, and hence also D, will be 
constants. We have merely R = Dr = Ds/h, and 


T (¢T 
De = ('e@o@)C@-adtda......., (5.1) 
Since w (x) = Oforx > T, we have further 
T 
Wit+r)-WO={ wOne+n-H —nt-Hlde=0 


fori > 7. Hence we have 

THeoreM II. If the originator is an arbitrary periodic function, 
the resultant will be periodic after equilibrium has been reached, and 
the length of its period will be equal to the length of the period of the 
originator. 

This condition, realized from the point of time 7 on, might be 
termed periodic equilibrium. For the constants h and H we have 


1 c+p 1 +p 
h=={ n(r)dr, H--{ W (r) dr 
P de P de 
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where c > T but otherwise arbitrary. Furthermore, from (4.2) 


T 
H=n{ w (€) dé 


whence 

THEOREM III. If the originator is periodic, the actual mean value 
of the resultant (after equilibrium) is equal to the constant value 
which the resultant would have had (after equilibrium) if the originator 
had been constant and equal to its actual mean value. 

For the function C (6) = C (4, 2), now independent of z, we have 


a C (6) = (ue +8) ry Oe... fo. aoe (5.2) 


where 
c+p 
a= \ u?(r) dr 


Since a = 0 is a trivial case, we may assume a > 0; C (4) is evidently 
independent of the constant c in (5.2), and by integrating (5.2) over 


5 we find 
c+p 
\ C (6) ds =0 


Furthermore, C (4) is an even periodic function with period p: 
C (6 + p) = C (8), C (6) = C (—6). Hence C is symmetric about 
6 = kp, where k is an arbitrary positive or negative integer or zero; 
C is even symmetric about 6 = kg, where p = 2q, for a C (¢q + 4) = 
a C (—q —8) is equal to 


c+p c+p 
i) u(r —q—a)u(s)dr = { u(r +q—8)u(r)dr =al (g —5) 


If u possesses derivatives of all orders, then (5.2) gives by v repeated 
integrations by parts, 


(m) es oye T? oe by ) 
aC (6) = ( yf U (r +6) u(r) dr.... (5.3) 


From this, we have form = 2p +1,» = p(p = 0,1,....), 
ic + 


2a Cert (0) oe (— yf ° alu” (r)]? =0 
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because [u'” (r)]? is periodic with period p. Hence C+»(0) 
= 0, since a >0. For m = 2p, »v = p, we getC® (0) = (—1)° a, 
where 


1 (c+? 
wm, = {ku Gide 
Qa c 
We thus obtain the Taylor series (assumed convergent) 
4 €3 PS eee 5.4 
(6) Sf ) 2p)! (5.4) 
Putting m = 1,v = 0, and m = 1,» = 1, in (5.3), and adding, we have 
c+p 
200’) =f [u’ (r +6) — wu’ (r — d) u(r) dr 


Hence C’ (gq) = 0. 
We shall next develop C in a Fourier series. Putting 


=> C (8) = C (2) 
toi u (t) = u (y) 


we have 
° ° c+2nr. c+2Q" : 
G(@) =c(-2),f (2) dz =0, f u(y) dy =0 


In the expansion 


C (2) =-5 Gm COS MX 


m=1 
assumed convergent, the coefficients a,, are given by 


c+2n . 
ran = { C (x) cos max dx 


Introducing the expression for C in terms of u, derived from (5.2), we 
get 


: c+2nr fc+2e ; 
ira = | i) u (x + y) u (y) cos mx dz dy... (5.5) 
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where 


c+2n 5 
a= fw way 


Assuming that « (y) can be developed in a Fourier series, let 6, and 
vu be the coefficients of this development, so that 


ti (y) - n(2) é =f nG)as 


= > (6, Cos wy + y, SIM wY,........... (5.6) 


Now substitute for u (x + y) in (5.4), leaving u (y) as it stands; 
expanding the terms cos (ux + uy) and sin (ux + wy), we see that in 
the integration with respect to x all terms containing sin u x. cos mx, as 
well as all terms containing cos yu x. cos mz, m + uy, will vanish, while 


+2 
f cos? mz dx = « 
¢c 


Introducing, next, the Fourier development of u (y) and integrating 
with respect to y, we get 


& A, = « (8), + 7,) 
From Parseval’s Theorem”, we have 
. aint: “3 “4 2 2 
a -{ u (y)dy=r = 6 +1,) 


Hence, finally 





C (6) = = dm 008 ( 


m=t 


=m) 


Pp 
where 


_ fn + Ym _ 

z (8) + ,) 

8, and y, being defined by (5.6). Since all the a,, are positive, it 
follows again that a; + a2 +.... = C (0) = 1 is the greatest value 


12 WHITTAKEK and Watson. Modern Analysis, 4ed., p. 182. 
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that C (6) can assume. Summing up our results concerning C (6), we 
have 

THEOREM [V. Let wu (¢) be an arbitrary periodic function with period 

ctp 
p = 2q, and such that | u(r) dr = 0. Then the function C (4) 
: ctp 
defined by (5.1) is also periodic with period p, and i) C (6) di = 0. 
If k is an arbitrary positive or negative integer, or zero, C (5) is sym- 
metric about all the points 6 = kg, and the first derivative C’ (6) 
vanishes at all these points; C (6) attains at the points 6 = kp its 
greatest value, viz., unity. The Taylor and the Fourier expan- 
sions of C (8) are given by (5.4) and (5.7), provided these series are 
convergent. 

This shows that even though the fluctuations of the periodic origina- 
tor n (t) be quite arbitrary, the function C (4) which occurs in the 
damping coefficient is subject to conditions that put definite restric- 
tions on the character of its variation. 

It may be seen that cos 2x 6/p is a function which satisfies the condi- 
tions to which C (6) is subject. 


6. THe DamPpiInG COEFFICIENT IN THE Periopic CAsE 


Substituting the Fourier expansion (5.7) in (4.8), we get 
© + 2 + @ 
D: = am {I f w (€) cos kn ede +[f w (é) sin mwede| b..(6.1) 


m=1 eJ—- @ 


where X,, = 2x m/p. For convenience, the integration has been 
extended over the interval (— ©, + ); this is allowable, since by 
definition w (x) = Oforxz < Oandz > T. 

We shall now introduce the expansion of w (x) in a series of Hermite 
polynomials or parabolic-cylinder functions: Let a and b be two 
arbitrary constants, and put z = x—a; then 


= k 
w (x) =w(z+a) =f(z) = ze Det (6.2) 


18 ARNE FisHEer. Mathematical Theory of Probabilities, vol. I. The convergence of 
this series is discussed in a recent paper by Cramir, On some classes of series used in 
mathematical statistics. Den Sjette Skandinaviske Matematikerkongres i Kébenhavn 
pp. 399-425. 1926. See also T. Kamepa, Theory of Generating Functions and its appli- 
cation to the theory of probability, Journ. Faculty of Science Imp. Univ. Tokyo, Sec. 1, 
1 (1): 1-62. 1925. Proc. Tokyo Math.-Phys. Soc., (2) 8: 262-295, 336-360. 1915. The 
Hp are tabulated by Jérgensen, Undersdgelser over frequensflader og korrelation, Kében- 
havn. 1916. Cf. Whittaker and Watson, Modern Analysis, chap. xvi. 























r: 


oct. 4,1928 MEIDELL: DAMPING EFFECTS AND EQUILIBRIUM 453 


in which 


1 2 
" (z) =¢ (z) = aren (- <<) 


f 
o.() = Se =(- 07H, (Z)o@ 


, 


d 
H a ce » _ 22/2 — 24/2 


hw of” f@ dH, (7) ere ithe (6.3) 


the H, being the Hermite polynomials. Substituting (6.2) in the 
general term of (6.1), and expanding cos \ (2 + a) andsin \ (z + a), 
we see that it is necessary to compute the two integrals 


+o 
Ep -{ cos \ 2° 9, (2) dz 


+o 
K, -{ sin \ z° y, (2) dz 


in which for brevity we have writtenXk=An. Since ¢g,(z) =(—1)’¢,(—2z), 
we have J2,,: = K2, = 0. Integrating by parts, we have from the 





second of (6.4), Ko,4:1 = —AJ2. To determine J2,, 2p integrations 
by parts of the first of (6.4) give 
3. ao ee ( =.) d 
ae ee _, C08 d 2 > exp 2 ha) 
d b)? 
= (— 1)’d”’ exp [- { » | 


hence 


» b)? 
Kiyar = (— DPt!N et exp | — 22] 


Using this, we see that the general term of the sum (6.1) will be 





Me gt el eee | st picky 
: J,,| ers ae = exp [- (Am 5)? [B,, + C,,I 


o=0 (2p)! 
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where 
8 k 2 7 
ba 2 (- 9 ots ( 7) 
dt aii — 
Pi k. 2a m\?r! 
ono 5 (a ayt ghee (Bem 
Te GT epi Y | 





The k, have the values given by (6.3); if desired, the k, can also be 
expressed in terms of the semi-invariants or the moments of w (2); 
the explicit expression for k, for any v in terms of semi-invariants has 
been given by Frisch." 

The expression for D is finally 


@ 


D? = = a, exp (— x? m®) - (BE +C*).........(6.6) 
where «x = 27 b/p, a» being given by (5.8), and B,, and C,, by (6.5). 
In numerical work, the arbitrary constants a and b would usually 


+o 
be disposed of by putting a equal to the mean M = \ éw(€)dE, 


+o 
and b equal to the stained deviation \¥? j= i (¢ — M)o dé. 


If this be done, k, and ky will vanish in (6:5). 
Since all the a,, in (6.6) are positive, and a, + a; +.... = 1, we have 


where Q is the maximum of 
Q,, = exp (— «° m’) - (Bi + C) 


for positive integral values of m. Hence we have the interesting 

TxeorEeM V. The limit (6.7) for D? involves the period p, but other- 
wise is completely independent of the originator. In particular, the 
magnitude of the fluctuations of the originator does not influence the 
limit of the damping coefficient if the originator is periodic with pe- 
riod p. 

If, for instance, the distribution function (i.e., the normalized 
distributor) is the normal, or Gaussian, error function, we have B,, = 1 


4 RaGNaR Friscu. Sur les semi-invariants et moments employés dans l'étude des 
distributions statistiques. Skrifter utgitt av Det Norske Videnskaps Akademi, 1926: 23. 
Oslo. 
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and C,, = 0; hence Q-= exp [—4(xc/p)*], where o is the standard 
deviation of the normal error function which represents the normalized 
distributor. Consequently 


2 
D < exp | - 2(=*) | Dy Suave s ewe (6.8) 
We therefore have 


TuHeoreM VI. If the originator is periodic with period p, and the 
distributor is a normal error function with standard deviation ¢, the 
damping coefficient is not greater than exp [—2 (x«/p)?]. 

In other words, the ratio of « and p hassuch a powerful influence on 
the magnitude of the fluctuations in the resultant that when the 
ratio ¢/p becomes small a practically perfect equilibrium—in the 
sense of a small value for D—will be brought about, regardless of 
the magnitude of the fluctuations in the originator. The damping 
effect is such that, e.g., if « = p, the value of D is less than 3.10-°. 
In such cases, it will rarely ever be necessary for practical purposes 
to take into account the fluctuations of the originator—the resultant 
may, under such circumstances, be computed as if the originator were 
constant and equal to its actual mean value. 

In view of Theorem V, the extremely low value of D found in the 
type of case just considered renders it reasonable to expect that a 
heavy damping will likewise take place when additional terms of the 
expansion (6.2) have to be taken into account. In each special case, 
the exact criteria have to be derived from (4.8) or (6.6) — (6.8); and 
the computed value of D has, of course, always to be viewed in relation 
to the nature of the particular problem under consideration. 

It may be expected that a closer study of the more general cases, 
when the underlying trend is not constant, would reveal, in substance, 
the same characteristic features in the damping effect. 


PALEOBOTAN Y—A palm fruit from the Miocene of western Panama.' 
Epwarp W. Berry, The Johns Hopkins University. 


The Tertiary flora of Central America is very incompletely known 
so that considerable interest attaches to any additions that tan be 
made to our knowledge of it. The palm fruit described below is such 
an addition. Unfortunately it has not been possible to assign this 


1 Received May 26, 1928. 
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fruit to an existing genus, but in extenuation it may be mentioned that 
not only are the recent palms of that region very incompletely known, 
but the majority are not represented by fruits in existing collections. 

Twenty-two genera of palms are enumerated in Standley’s Flora of 
the Canal Zone,? but all of these except the genera Acanthorhiza and 
Geonoma are quite different from the fossil, and there must be many 
species and even genera in Central America as a whole which remain 
undiscovered. 

The present fossil is referred to the form genus Palmocarpon, and 
may be described as follows: 


Palmocarpon geonomoides Berry, n. sp. (see fig. 1) 


Fruit small, symmetrical, nut-like, nearly spherical—being 18.5 millimeters 
long and 17 millimeters in diameter. Lacking a raphe or any trace of a 
tripartite division, or of the basal pores of the 
cocoid palms. Fibres simple and flat, parallel with 
one another to the extreme base, not over 1.25 milli- 
meters in maximum width, continuous from the base 
to within a millimeter or two of the apex, where their 
acute points are slightly raised as if free tipped 
in life; united by their edges to form a pseudo-shell. 


The single specimen is a sandstone cast with 
, a very slight film of carbonaceous matter in 

Fig. 1.—Palmocarpon places. The general form is simulated in ex- 
geonomoides Berry, n.8p., isting species of Oenocarpus, Geonoma, Acan- 
Miocene of western Pan- 2 ‘ 
Praia thorhiza, etc., and the fibres suggest comparisons 

with existing species of Rhopalostylis and the 
outer coat of Astrocaryum. The fibres are, however, more regular 
than in Astrocaryum, especially toward the base, pores are lacking, 
and there is no trace of an inner fibrous layer which could hardly 
have disappeared during fossilization. Nor does it seem possible 
that the surface in life could have had several layers of imbricated 
lanceolate flat fibres without some trace of them having remained on 
the fossil. The fact that the specimen is fully inflated and perfectly 
round in transverse profile shows that it did not grow in a crowded 
inflorescence. 

Because of the uncertainties of a positive generic identification and 
the lack of adequate recent comparative material it is referred to the 
form genus Palmocarpon instituted by Lesquereux in 1878 for palm 
fruits of uncertain generic affinities, the specific name being indicative 





* Stanpiey, P.C. Contr. U.S. Nat. Herb. 27: 93-100. 1928. 
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of my opinion that it belongs to the tribe Geonomeae. The specimen 
comes from half a mile northwest of Zembrano, Province of Chiriqui, 
in western Panama, and was collected by James Terry Duce, and 
presented to the National Museum by T. D. A. Cockerell (U. 8. N. M. 
no. 37193). It is labelled as having come from the Uscari formation, 
which on the Caribbean side of the Isthmus lies beneath the Gatun 
formation, but which has not been definitely recognized in Chiriqui 
Province. The age is undoubtedly Miocene. 

Silicified palm wood was described by the writer® in 1918 from the 
Cucuracha formation (Oligocene) in the Gaillard Cut; and an Iriartea- 
like fruit was described‘ from the Gatun formation one and a half 
miles northeast of Gatun in 1921. These with the foregoing comprise 
all that is at present known of the Tertiary palms of Panama. No 
palms were represented in the small Tertiary florule described from 
Costa Rica in 1921,° but palm rays are not uncommon in collections 
from the Isthmus of Tehuantepec in southern Mexico.‘ 


PROCEEDINGS OF THE ACADEMY AND AFFILIATED 
SOCIETIES 


PHILOSOPHICAL SOCIETY 
972D MEETING 


The 972d meeting was held at the Cosmos Club March 17, 1928. 

Program: F. L. Mouuer: Recombination of atomic ions and electrons. 
Spectra result from recombination of atomic ions and electrons. When a 
thermionic discharge tube designed to favor high concentration of ions was 
operated in caesium vapor with relatively high current and low voltage, 
continuous bands appeared at the limits of each series, in general agreement 
with theoretical predictions. Electrical measurements of discharge condi- 
tions have been made by the Langmuir probe wire method. The current 
voltage curves to the probes show that the electron velocities have a strictly 
random or temperature distribution with very low average speeds ranging 
between 0.2 and 0.3 volts. The ion concentration is of the order of 10" per 
ec. On the basis of data on the intensity distribution beyond a limit and the 
electron velocity distribution one can compute the chance of recombination of 
an electron into the corresponding energy level as a function of the electron 
speed. This probability decreases very rapidly as the speed increases. 
From this can also be derived the probability of absorption of an atom in this 
state. The continuous emission bands beyond the 2P and 3D limits of 
caesium give similar absorption curves which drop rapidly from the limit to 


* Berry, EowarpvW. U.S. Nat. Mus. Bull. 108: 24. 1918. 
‘Berry, EpwarpW. Proc. U.S. Nat. Mus. 69: 21-22. 1921, 
5’ Berry, Epwarp W. Idem. 169-185. 

* Berry, Epwarp W. Idem. 62 (Art. 19): 5, pl. 5, fig. 4. 1923. 
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higher frequency. The form of curve is apparently quite different from the 
theoretical, which decreases approximately as the 4th power of the wave 
length. It resembles the curve for the 1S state derived by Foote and Mohler 
from photoionization measurements. (Author’s abstract.) 

R. M. Lancer: Dispersion and quantum theory. The interaction between 
radiation and matter can be expressed in terms of a property of the atom, 
namely the electric moment. The expressions for the rate of emission and 
absorption of light were given and discussed. The connection with the 
duration of emission or life of an atom in the excited state was mentioned. 
The relation between electric moment and absorption coefficients as well as 
refractive index was pointed out. Slides were shown of experimental results 
which proved that the positions of anomalous dispersion, absorption and 
emission were identical. This fact was decidedly in disagreement with the 
picture of the atom developed by Bohr and the discrepancy was one of the 
chief sources of dissatisfaction with this older form of quantum theory. The 
new quantum mechanics gives rise to a picture of the atom in which there is a 
movement of charge which matches in frequency the light given out, or 
absorbed. The positions of anomalous dispersion are then those to be 
expected from general electromagnetic theory and coincide as found in 
experiment with the absorption lines. We can therefore have more faith in 
the new model of the atom given by Schrédinger and it becomes interesting 
to calculate some of its properties. The distribution of charge in normal 
hydrogen was computed and a section of it was portrayed as the concentration 
of white dots on a black background. The atom looked like a spherical 
globule of charge with density decreasing gradually from the center outwards. 
The sphere with radius equal to the radius of the atom on the Bohr picture 
contained only half the charge. The considerable concentration outside this 
radius gives the impression of stray fields and fits in with our notions of the 
chemical notions of hydrogen. The excited state of hydrogen 2S was shown 
on a similar slide. It is also spherically symmetrical but the charge is not so 
concentrated, and is distributed as a diffuse shell concentric with the inner 
globule. In general S states of quantum number n would have n shells. 
The lithium atom is similar to 2S hydrogen, the sodium atom to 3S hydrogen, 
and so on. Normal helium was also shown and contrasted with hydrogen. 
The density at the center is very much higher but diminishes much more 
rapidly with radius so that practically all of the double electronic charge is 
within 10-* cm. of the nucleus. This means there is very little stray field as 
we should expect from the fact that helium is an inert gas. The helium 
structure is characteristic of the outer shell of the elements of the second 
column of the periodic table, e.g., Hg and Cd. Another property which can 
be calculated with the help of the quantum mechanics is the atomic absorption 
coefficient. The total absorption of many molecules will involve an integra- 
tion and not a simple multiplication by the number of molecules. 

The mean life can also be computed and for the atom in the second quantum 
state has a value of about 1.58 XK 10-*seconds. (Author’s abstract.) 


973D MEETING 


The 973d meeting was held at the Cosmos Club March 31, 1928. 
Program: C. J. Davisson: Reflection and diffraction of electrons by a 
crystal of nickel. It has been found in recent experiments by the writer and 
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Dr. L. H. Germer that a stream of electrons is regularly reflected from a face 
of a nickel crystal. This phenomenon could be understood if electrons were 
waves, but is incomprehensible in terms of atoms and electrons and their 
interactions as we have pictured them. Reflection is not observed from a 
polycrystalline target, which suggests that if electrons are waves, their 
wave-lengths, like those of x-rays, must be comparable with the distances 
between atoms in solids. The reflection from the crystal is selective in speed 
of bombardment; intensity maxima occur at nearly equal intervals in this 
variable. This would be comprehensible if electrons were waves of wave- 
length inversely proportional to their speed. These results suggest by 
analogy with x-ray phenomena, that electrons will be regularly reflected also 
from sets of atom planes which do not lie parallel to the surface of the crystal. 
Beams of electrons of this type are found issuing from the crystal at 
critical speeds of bombardment, but their directions are not those of regular 
reflection from the principal atom planes. The difference between x-ray and 
electron diffraction in this respect is due apparently to a difference in the 
refractivity of nickel for the two kinds of radiation. It is shown that if this 
is the case the wave-length of the diffraction beam will nevertheless satisfy 
the plane grating formula with respect to the atomic plane grating lying 
parallel to the crystal face—files of atoms serving as the lines of a grating. 
This formula which does not involve the refractive index of the crystal is used 
to calculate electron wave-lengths, and all values so found agree, within the 
limits of accuracy of the measurements, with the corresponding values of 
h/mv—Planck’s constant of action divided by the momentum of the incident 
electron. They agree, that is, with the wave-length associated with a freely 
moving particle in the undulatory mechanics introduced by L. de Broglie. 
The regularly reflected beams cannot be used for calculating wave-lengths. 
The Bragg formula cannot be used because the refractive index of the crystal 
is not unity, and the plane grating formula cannot be used because the beams 
are of order zero. Assuming h/mv to be the wave-length of the incident 
beam the data of regular reflection may, however, be used to calculate indices 
of refraction. When this is done it is found that the index of refraction of 
nickel for electrons decreases from about 1.15 to 1.01 as the bombarding 
potential is increased from 60 to 600 volts. (Author’s abstract.) 

H. L. Curtis in an informal communication called attention to the close 
agreement between the recently measured velocity of light and value of the 
ohm, and the calculations of Rosa and Dorsey twenty-one years ago based 
upon their measurements of the ratio between the electromagnetic and electro- 
static units of electricity. 


974TH MEETING 


The 974th meeting was held at the Cosmos Club April 14, 1928. 
Program: L. W. Kepuart: Plant hunting through East Africa. 
H. E. Merwin, Recording Secretary 
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SCIENTIFIC NOTES AND NEWS 


In connection with the Thirteenth Annual Meeting of the Optical Society 
of America an exhibition of optical instruments, apparatus and products will 
be held under the joint auspices of the Bureau of Standards and the Optical 
Society, from October 31 to November 3. 


Dr. A. 8. Hrrcencock has returned from Newfoundland and Labrador, 
where he has been studying and collecting grasses. In Newfoundland he 
collected at Port-aux-Basques, St. Georges, Corner Brook, Little Harbor, 
Grand Falls, and St. Johns, and in Labrador at Battle Harbor and Cartwright. 


Dr. C. G. Assot, Secretary of the Smithsonian Institution, returned to 
Washington on September 23 after a successful expedition to Mount Wilson, 
California. Bolometric observations of the infra-red solar spectra and radio- 
metric observations of the spectra of many stars formed his principal work. 
The results, which have not yet been fully worked up, appear to be very 


gratifying. 


Mr. M. W. Stirina, Chief of the Bureau of American Ethnology, attended 
the meeting in New York City of the XXIII Session of the International 
Congress of Americanists, which convened on Monday, September 17. Mr. 
Stirling was appointed by the State Department to represent the United 
States Government at the Congress. He returned to Washington Septem- 
ber 24. 








